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PACS 81 . 05 . ue - Graphene 

PACS 72 . 80 . Vp - Electronic transport in graphene 

Abstract. - We study the effect of extended charge defects in electronic transport properties of 
graphene. Extended defects are ubiquitous in chemically and epitaxially grown graphene samples 
due to internal strains associated with the lattice mismatch. We show that at low energies these 
defects interact quite strongly with the 2D Dirac fermions and have an important effect in the 
DC-conductivity of these materials. 



Introduction. — Graphene crystals isolated by the 
exfoliation method [I] are high quality films, with high mo- 
bilities even on a SiC>2 substrate. While exfoliation works 
well for the study of the fundamental physical properties 
of graphene [2] it is not a scalable process useful for tech- 
nological applications [3j. At the present time, the most 
promising scalable growth methods of graphene films are 
either based on epitaxial growth on SiC 4 or on chemical 
vapor deposition (CVD) of graphene on metal surfaces [5j. 

Graphene growth on crystal substrates, independently 
of the process, is subject to strain due to lattice mismatch 
between graphene and the substrate. The strong sigma 
bonding between the carbon atoms makes the graphene 
lattice very stiff (the spring constant is of order of 50 
eV/A 2 ) and hence in-plane deformations are energetically 
costly. Because graphene is the ultimate example of a 2D 
film, the strain can be release by two main mechanisms, 
namely, either by going out of the plane or by the recon- 
struction of the chemical bonds. 

By exploring the third dimension the graphene film pays 
the energetic price of bending (the bending rigidity is of 
order of 1 eV) and the loss of interaction energy with the 
substrate (which is usually a mix of covalent bonding and 
van der Waals interaction). In certain cases, the energetic 



price of forming wrinkles and blisters in the graphene film 
is smaller than the price of creating structural defects such 
as pentagons, heptagons, and octagons. STM studies have 
shown that epitaxial graphene grown on 6H-SiC actually 
bends and buckles as a result of the compressive strain [6] . 
The same effect is observed in samples grown by CVD on 
Cu [7]. In other cases, when the the interaction between 
graphene and the substrate is strong, it is energetically 
preferable to reconstruct the bonds with the formation 
of lines of defects. This is what happens, for instance, in 
graphene grown by CVD on a Ni surface [8] . The surface of 
the film reveals the presence of extended one-dimensional 
defects. These defects are lines of periodic cells made of 
two pentagons followed by one octagon. 

The general conclusion is that the intrinsic 2D nature 
of graphene makes the presence of one-dimensional ex- 
tended defects in artificially grown graphene samples a 
norm. The formation of extended structural defects in 
graphene has strong consequences for the electronic prop- 
erties. On the one hand, localized bending and strain 
can lead to the appearance of strong "pseudo-magnetic" 
fields that create localized Landau levels. This is the main 
idea behind the concept of strain engineering [9]. On the 
other hand pentagons, heptagons, and octagons, act as 
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Fig. 1: Tapping mode atomic force microscope (AFM) image 
of a CVD graphene film transferred on Si02- Different ther- 
mal expansion of the Cu foil and the graphene sheet result in 
the formation of a few nm high ripples. Locally cracks can 
form during the transfer process and occasionally one is left 
with PMMA residues. Inset: Hall bar device used to measure 
the conductivity of CVD grown graphene (scale bar is 5 /mi). 
Representative data is given in Figs. [2] and [5] 



donors and/or acceptors relative to the normal benzene 
ring configuration, and hence also leads to charge local- 
ization [To] . All these mechanisms are examples of the 
more general concept of self-doping in sp 2 bonded carbon 



11 . The existence of the localized states at the defect 
line lead to transfer of charge from the bulk of graphene 
to the defect — charge accumulates at the extended de- 
fect, creating charged lines. Such type of scattering lines 
can be the limiting scattering mechanism of the electronic 
mobility in these graphene films. Furthermore, experimen- 
tal studies show that these extended defect lines intercept 
each other at random angles, forming irregular polygons 
with edges showing stochastic distribution of length. A 
study of the electronic scattering in chemically produced 
graphene needs to take into account such types of random 
distributions in order to correctly describe the effect of the 
extended charge defects in DC-transport. 

This Letter is organized as follows; we begin by outlin- 
ing the experimental procedure used to produce the CVD 
graphene samples and briefly characterize the extended 
defects seen in the microscopy studies. Then, the main re- 
sult is presented and tested against DC-conductivity data 
for several CVD graphene samples. In the remainder we 
discuss in detail the model of extended charge defects, the 
effective scattering potential and its contribution to the 
semi-classical conductivity. Last, finer points of our cal- 
culation, such as the effect of disorder in the length distri- 
bution of the defects, and the electron-hole asymmetry of 
electronic cross-sections due these defects are discussed in 
separate sections. 

Outline. — The theory described below has been used 
to interpret the transport data of graphene grown using 
the roll-to-roll method fl2|. An AFM image of a CVD 
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Fig. 2: Conductivity measured in CVD synthesized graphene 
at T=3.5K (data points shown in blue) [H] is fitted to Eq. [I] as 
function of V = V g — Vmin (dashed line). The optimal parame- 
ters are found to be W ~ 14.6nm; 7 ~ 2.14 X 10 u cm~ 2 (V > 0) 
and 7 ~ 2.28 x 10 11 cm~ 2 (F < 0). The inclusion of midgap 
states (i.e. resonant scatterers) is shown in thick lines [red 
(V > 0) and green (V < 0)] and modifies 7 to 1.15 x lf/V -2 
(V > 0) and 1.27 x 10 n m" 2 (V < 0). The midgap parameters 
are n s = 1.6 x 10 11 cm -2 and R = 2ao, with ao = 0.14nm; the 
experimental data was shifted as to have a minimum of zero 
at the Dirac point (Vmin —5.5V and a m i n ~0.1e 2 //i). Fits to 
other CVD samples are shown in Fig. [5] 



tended line defects, few nanometers long, can be seen. 
Sample preparation and measurement were performed us- 
ing standard methods: graphene samples are grown on Cu 
substrates by chemical vapor deposition (CVD) 12 . To 
characterize CVD graphene samples, standard Hall bars 
are patterned by e-beam lithography, followed by thermal 
evaporation of Cr/Au (5/25 nm). An additional e-beam 
lithography step followed by O2 plasma etching is per- 
formed to define Hall bar device as shown in the inset of 
Fig. [T] Measurements are performed as a function of tem- 
perature down to 3.5 K using standard lock-in amplifier 
techniques. 

The central result of this work is an expression for 
the semi-classical DC-conductivity of graphene due to ex- 
tended charge defects, 
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graphene film transferred to Si02 is shown in Fig. [T] ex- 



where rii is a finite density of effective extended de- 
fects (see next section), made of lines with charge q% and 
G~ l {x) = J dO (1 - cos 6) (tt - Of [5 + 4 cos (x-x cos 6)]. 
For sake of simplicity, in deriving the latter equation, we 
have considered a graphene structure constant a = 1/2. 
The Fermi momentum Uf relates to the electronic car- 
rier density n c according to fcj? = ^/7m c , which can be 
controlled by the application of a back-gate voltage, V g , 
after transferring the graphene sheet to a dielectric sub- 
strate, typically silicon oxide, 300nm thickness, for which 
one has n c — 7.2 x 10 14 V g (SI units). To test our results, 
we used experimental data of conductivity measurements 
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according to: 



Fig. 3: Left — Schematic of a small area of CVD graphene with 
charged lines with several length; the area (in yellow) encom- 
passed by 3 lines defines the intersection zone of the effective 
extended charged defect. Right — The effective extended de- 
fect is built, from three lines initially at x = 0, translating one 
line (e.g. £3) along the x-axis by W and rotating the lines by 
ipi about the 2-axis. In the picture (pi — to make it easy to 
visualize; the origin of the coordinate system is represented as 
a black dot. 



in several graphene samples grown by CVD on Cu, using 
the roll-to-roll method 12 . The fitting parameters in our 

niqf/e 2 . The former has units 



theory are W and 7 
of length and is roughly the mesh size originated by the 
intersection of defect lines, whereas the latter measures 
the scattering strength of the extended defects. Figure [2] 
shows a high-quality fit to the CVD data. For moderate 
to high carrier densities, we find that Eq. [T] fits perfectly 
well the data, indicating that linear charged lines formed 
during the CVD growth process can play a key role in 
DC-transport at finite electronic densities. The line sep- 
aration was found to lie within the range 10-15 nm in all 
samples. We also found that the inclusion of strong (res- 
onant) short-range disorder improves the agreement with 
the experimental data for lower carrier densities, more pre- 
cisely, \V g — V m i n \ < 10V. (Recall that resonant scattering 
has recently emerged as one of the main mechanisms limit- 
ing DC-transport of non-suspended graphene [13[|14] .) In 
what follows, we give a detailed description of our model 
of extended charge defects in graphene. 



Theory. — We start by characterizing the scattering 
potential created by a single extended charge defect. We 
take the charged defect to be an infinite line along the y 
axis, with linear charge density A/, embedded in graphene, 
which is lying in the xy— plane. Its 3D charge density 
has the form pu nc (r) = \i5(x)5(z). Basic electrostatics 
predicts that single line of charge in vacuum produces a 
logarithmic potential in space. Clearly, embedded in a 
metal, the potential is modified by screening effects. The 
screening can be taken into account within the Thomas- 
Fermi (TF) approximation [15] . The TF assumes that the 
local electronic charge density, p(r), changes due to the 
effective potential, </> e ff(r), created by the extended defect 



6p(r) ~ -ep(E F )e(j> cS (r), 



(2) 



where p{Ep) = 2kp/ {%hvp) is the bare density of states 
per unit of area (spin and valley degeneracies included) 
and kp (vp) is the Fermi momentum (velocity). 

The effective potential, (f> e g, is determined by Poisson's 
equation: V 2 </> cff = - (p linc + Sp) /e d , that is, 
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where e d is the dielectric constant of the medium [notice 
that Sp = — e 2 p(Ep)(j) e ff(x, z = 0)]. This equation can be 
solved by Fourier transform followed by an integration over 
the z coordinate. The form of the potential in momentum 
space is 

&ff(gs) — ;r~ III ■ ( 4 ) 
££d \1x\-r QTF 

where qTF = ^akp is the TF wave vector, with the 
effective graphene's structure function given by a = 
e 2 /{A-Ke d hvp ). We note that (f> e &(x) shows a logarith- 
mic divergence at the origin and asymptotic behavior 
(qrFX » 1) given by 



4>es{x) -> hi [2e d -K{q TF xy 



(5) 



In contrast with the screened Coulomb potential created 
by a point charged impurity, the potential of a charged 
defect line does not decay exponentially away from the 
scattering center, and hence can lead to a strong effect on 
electronic transport, as we show below. 

Atomic force microscopy of small areas of graphene 
(^0.1/im 2 ) shows extended defects with several lengths 
and orientations [5j. We model these extended defects as 
straight lines that intersect forming polygons (see Fig. [3]) . 
The electronic scattering is determined by the polygons 
formed by the defect lines (see later). The network of 
defects is built from a number of base lines (labeled li), 
initially lying along the y-axis, in two steps: translating 
each line U by a vector Rj = (xj,t/j), and finally rotating 
them about the z-axis by random angles {ifi}. Through- 
out the paper we assume that the size L of these lines is 
much larger than other length scales in the problem. For 
N base lines the screened potential reads: 
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where Q{<Pi) = q x cos ipi — q y smipi is the projection of 
the wave vector q onto the direction perpendicular to the 
line li, defining the direction of momentum transfer in 
an electron-defect scattering event, and qi = XiLi is the 
charge of line li. [To obtain this result, we approximated 
the 2D electrostatic potential of a single line (e.g. oriented 
along the y axis) using </> c ff(q) = / dycf> cfi (q x )9(L/2 -\y\), 
where 9{y) is the Heaviside function and (j) e g(x) is the 
potential of an infinite charged line embedded in graphene 
(Eq-i).] 
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Model of the extended defect. Transport through 
nano-electronic graphene devices with extended charge de- 
fects is tackled here from the point of view of a single 
effective defect "cell", from which a prototype network 
of charged lines can be reproduced. The effective ex- 
tended defect is made of three lines, with equal lengths, 
which intersect forming a triangle (see Fig. [3]); choosing 
R-i(2) = (0,0) and R3 = (W, 0), introduces a new length 
scale, W, the amount of translation of I3, here referred 
to as line separation. Its scattering potential (denoted by 
0A,cff) is given by Eq. (|6| setting N — 3. Figure G] shows 
how many such defects give rise to all sort of polygons. 

The large-distance behavior of the potential due to a 
charged line in graphene (see Eq.[4]) renders the Born series 
particularly suitable to compute scattering amplitudes. In 
the first Born approximation (FBA), the elastic scattering 
amplitude for massless Dirac fermions in 2D reads 16 



f(0) = 




A, off 



(q)>, 



(7) 



where p (|p| = kp) is the wave vector of the incident 
electron (we choose p = |p|e x ), q = p' — p the trans- 
ferred wave vector (p' stands for the "out" wave vector), 
6 = Z(p',p) is the scattering angle and (0A,eff(q)) is the 
scattering potential conveniently averaged as to include ro- 
tational disorder (disorder in W will also be considered); 
the factor = 1 + e lB originates from graphene's Berry 
phase precluding electrons from back-scatter [17]. The 
main result ([TJ) is obtained within the Boltzmann approach 
using the FBA, 

at = {\/2)e 2 v 2 F p{E F ) Tu (8) 
where the relaxation time is given by I/77 — nivp J d6(l — 
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Rotational disorder. Within the model depicted in 
Fig. |3j rotational disorder can be taken into account by 
averaging the scattering potential </>A,eff over the angles 
{ifii}. The scattering potential due to a single A extended 
charge defect was averaged using a uniform random dis- 
tribution of angles: 



(0A,eff(q)) = 7T~ 



JJ<fy>i^A,eff(q)- 



(9) 



The influence of the particular lines orientation on trans- 
port depends strongly on the screening length. For 
graphene on top of SiC>2, the graphene structure constant 
a is expected to be around 0.5 and thus qxF ~ 2fcp. In this 
case, the variation of the extended charge defect potential 
(4>A,es) with the angles {tfii} is hindered by a TF wave- 
vector qtf with the same order of magnitude than Q(cpi), 
as direct inspection of Eq. Q shows. As a consequence, 07 
becomes little sensitive to the particular random distribu- 
tion adopted. At the Fermi level, the averaged scattering 
potential reads (see Appendix), 
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Fig. 4: An effective extended defect with triangular intersection 
(three lines in black) is replicated and translated according to 
the green arrow, producing a mesh containing several types of 
irregular polygons. 



Interestingly, the line separation W adds an oscillating fac- 



tor to the scattering amplitude [last term in ( 10 )], with the 



consequence that the familiar V-shape in the conductiv- 
ity as function of the gate voltage [l3j[l4] (which requires 
f(9) ru 1 / y/kp) will not be observed in samples with a high 
density of extended charge defects. This manifests into a 
change of curvature in a a vs V g plot, a bona fide signa- 
ture of scattering due to extended charge defects. This 
oscillating factor is essential to the high quality fit of the 
data, as shown in Figs. [2] and [5] 

Disorder in the mesh size. Experimental studies in 
CVD graphene show extended defects with many geome- 
tries, the fact that the A extended defect fits well the data 
indicates that such defect constitutes the dominant type 
of scatterer within our model. As discussed above, albeit 
for a ~ 1 /2 the conductivity is barely affected by the spe- 
cific orientation of the lines, it is very sensitive to changes 
in W, since this parameter measures roughly the defect 
mesh size, and hence is directly related to the scattering 
strength experienced by the electronic carriers. A care- 
ful inspection shows that a change of 5% in W is enough 
to deteriorate the fits. This parameter sets the departure 
from the a ~ n c (or equivalently, a ~ V) behavior, similar 
to that originated from charged impurities located in the 
substrate, to a more involved carrier-density dependence, 
namely that of Eq. ([IJ . 

The effect of disorder in W can be estimated by assum- 
ing a normal distribution p(W) with average line separa- 
tion W and variance dW 2 . The potential accounting this 
kind of disorder is obtained by replacing the exponential 



factor in Eq. 10 according to 



(11) 



where q x — kp (cos — 1) is the transferred momentum 
along the x axis; a finite variance dW 2 changes the 
previous results (i.e. without disorder in W) whenever 
(q x dW) 2 > 0(1). We find W/10 > dW > W/100, in 
all samples, showing that our model predicts that promi- 
nent extended defects form intersection edges with lengths 
centered around W ~ O(10nm) with small variance. 
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Fig. 5: Conductivity measured in CVD synthesized graphene 
at T=3.5K for two samples with different mobilities (blue dots) 
and respective fits to the semi-classical calculation with res- 
onant scatterers contribution included. Sample with higher 
(lower) mobility: W ~ f0.0(ff.5)nm; 7 ~ 5.38(4.84) x 



'(V < 0). The 
2 and R 



10 lu cm _/ (l/ > 0) and 7 ~ 4.63(5.56) x 10 lu cnT 
midgap parameters are n s = 1.5(3.0) x 10 cm' " ana ti = ao 
with ao = 0.14nm; the experimental data was shifted as to 
have a minimum of zero at the Dirac point: V m m — 8.0(5. 5)V 
and a min ~0.2 (0.1)e 2 /h. 



Other scattering mechanisms. So far we have analyzed 
the effect of extended charge defects in DC-transport in 
graphene. Notwithstanding, other mechanisms are in play 
which can provide important corrections to our model, 
especially in the regime of low carrier density, where the 
fit to Eq. [I] is less accurate (see dashed line in Fig. [2]) . 

We focus our attention on midgap states, presumably 
the most important scattering mechanism in mechani- 
cally cleaved non-suspended graphene samples at finite 
densities and not too high temperatures 13 Mj 18 



|ig. 

here also playing an important role as we will briefly see. 
Midgap states emerge due to resonant scatterers (RS), 
whose physical realization could be vacancies or adsorbed 
hydrocarbons in the surface of graphene. Indeed, we as- 
sume a typical value for the density of resonant scatterers, 
n s ~ 10 1:L cm _2 , and take the radius of the scattering disk 
R to be of the order of carbon-carbon distance. The cor- 
rection to the conductivity (Eq. [I]) is then calculated using 
Matthiessen's rule, er -1 = a^ 1 + erj 1 , where, 



2e 2 

k 2 F In 2 (k F R) 



hn 2 



(12) 



is the conductivity due to resonant scatterers 13 14 20 



The new fits are obtained by keeping W fixed from its 
previous value (i.e. with just extended charge defects con- 
sidered) and varying 7 — Fig. [2] shows that midgap states 
yield an important correction in the low to moderate den- 
sity regime. 

Electron-hole asymmetry. We finally discuss the ori- 
gin of the electron-hole asymmetry in DC-transport, high- 
lighted in Fig. [2] by using different colors to represent p 
and n conductivity. The asymmetry between holes and 



electrons mobility 



H = (l/e)da/dn c , 



(13) 



is clearly shown in most transport studies in graphene; in 
CVD samples \fx p — (inl/Hn is about a few percent. This 
asymmetry may have two distinct physical origins: (1) 
a potentially significant charge transfer from the metallic 



contacts to graphene 21 , and (2) scattering cross-sections 
sensible to the carriers polarity. 

In order to estimate the contribution of the extended 
defects to the conductivity asymmetry, we compute the 
next term in the Born series for a single charged line. The 
2 nd Born correction to the scattering amplitude Sf(9) is, 



Sf(0) 



2p T 



evph 



d 2 k - 



(p'-k)Gz>(k)<Mk-p) 



(14) 



/(*) '4ff(q)J (2tt) 
where the 2D Dirac propagator reads 

Gc(k) = v 2 F h 2 (|p| 2 - k 2 

and the screened potential (f> c g is given by Eq. [4] with 
XiL — qi. The differential cross-section \f + 5f\ 2 has now 
a term proportional to 4>%q, and thus is no longer invariant 
under a change of electrical charge sign e — > — e. For small 
density of charge \qi/e\ < 0.1, the conductivity in the 2 nd 
Born approximation, a^ 2 \ relates to the FBA value, <j^\ 
according to 



A 2 ) 



1 - g(a) 



(15) 



where g{a) reflects the importance of the dielectric 
medium. For graphene on top of silicon oxide we obtain 
5(1/2) ~ 1.2, entailing a negligible difference between the 
conductivity for carriers with opposite polarities as long 
as qi -C e. In general, for extended defects made of a 
single line of charge, the 2 nd Born approximation yields 
a global multiplicative factor in the conductivity; notice 



that the term inside brackets in Eq. (151 just depends on 



the charge of the defect and on the sign of carriers, qi/e, 
and therefore is constant throughout the entire range of 
carrier density. Indeed, one can absorb the correction from 
the 2 nd Born approximation into 7 — niqi/e 2 according to 

7 -> 7/[l -ff(a)<fi/e]. 

It would be desirable to perform the same calculation 

for a A defect. Unfortunately, however, in this case, the 
2 nd Born correction becomes much involved. Despite that, 
we can get some intuition by studying a specific configura- 
tion; we have performed a numerical evaluation of g(a) for 
aJV = 3 extended defect with </ 5 i(2)(3) = 0, R-i(2) = (0, 0), 
R3 = (W,0) (see Figure [3] for definitions of tpi and Rj) 
and W =10nm; a strong dependence of g(a) with kp was 
observed, more precisely, a variation of ~30% by increas- 
ing the gate voltage from IV to 50V. We leave as an open 
question whether the inclusion of the 2 nd Born amplitude 
in the calculation of the conductivity leads to a qualita- 
tive improvement of the fits of experimental data to our 
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model — this would elucidate about the precise amount of 
electron-asymmetry possibly induced by charged extended 
in CVD graphene. 

The experimental data shown in Figs. [2] and [5] disclose a 
small but perceivable change of mobility from p to n car- 
riers. In the light of the latter results, this suggests that 
the 2 nd Born correction cannot be too large in these sam- 
ples (discarding the remote possibility that a large second 
order correction is balanced with a charge transfer from 
the Hall probes). 

Conclusions. — Extended defects are prevalent in 
CVD graphene and arise in SEM studies as cracks with 
several sizes and oriented at random. Due to the self- 
doping mechanism 
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such cracks will act as charged 
scattering centers. In this Letter, we have studied the 
impact of such defects in the DC-transport properties of 
graphene films. By constructing a simple model of ex- 
tended defects, a semi-classical computation of the DC- 
conductivity was carried out taking into account disor- 
der on the extended defects geometry and the screening 
by graphene electrons. We have shown that charged ex- 
tended defects lead to a very distinctive dependence of 
DC-conductivity with carrier density compared to previ- 
ously studied mechanisms 13 14 . Our findings show that 



extended charge defects can play an important role in DC- 
transport of chemically grown graphene samples with a 
large density of cracks. 

Growing graphene via CVD is a very promising route to- 
wards scalable fabrication of two-dimensional high-quality 
carbon films. The understanding of the scattering mecha- 
nisms in chemically grown graphene samples is thus of up- 
permost importance to increase electronic mobilities cur- 
rently limited to ~4000 cm 2 -V -1 -s -1 . Given the stringent 
constraints on electronic mobilities required for techno- 
logical applications of graphene, our results show that the 
control of such defects can be of fundamental importance 
for further development of a graphene-based electronics. 
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Appendix. — The Thomas-Fermi renormalized elec- 
trostatic potential for the extended charge defect A de- 
picted in Fig. [3] reads, 



^A,eff(q) = 



qi_ \ 1 

2e<i [ \ <}x cos (fi qy sin ip\ | + (Jtf 
1 

+ 1 : 1 

\q x cos v?2 - q y sin ^ 2 1 + <7tf 



\q x COS ip 3 - q y Sin ip 3 1 + q TF 



(16) 



The potential depends on the relative orientation of the 
lines through the angles {fi}. To get a meaningful re- 
sult we have to consider a statistical distribution of such 
angles. For a > 1/2 the denominators in (16) are dom- 



inated by qxF and hence we can safely integrate the 
{<fi} dependences using a uniform distribution. Choos- 
ing q = kp(cos9 — 1, sin#), we obtain, 
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(17) 



with w a = y/2{8a 2 + cos0- 1), f° = |sin6»|, f* 
G±/(sin 2 §) 



5Q 

where 



±1 



G± = [cos (9- It 4a] 



±2a|sin0|. (18) 



Setting a = 1/2 in (17) one obtains Eq. 10 
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